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ABSTRACT 

Hot Jupiters, with atmospheric temperatures T > 1000 K, have residual thermal ionization levels sufficient 
for the interaction of the ions with the planetary magnetic field to result in a sizable magnetic drag on the 
(neutral) atmospheric winds. We evaluate the magnitude of magnetic drag in a representative three-dimensional 
atmospheric model of the hot Jupiter HD 209458b and find that it is a plausible mechanism to limit wind 
speeds in this class of atmospheres. Magnetic drag has a strong geometrical dependence, both meridionally 
and from the day to the night side (in the upper atmosphere), which could have interesting consequences for 
the atmospheric flow pattern. By extension, close-in eccentric planets with transiently heated atmospheres 
will experience time-variable levels of magnetic drag. A robust treatment of magnetic drag in circulation 
models for hot atmospheres may require iterated solutions to the magnetic induction and Saha equations as the 
hydrodynamical flow is evolved. 

Subject headings: 



1. INTRODUCTION 

Hot Jupiters are close-in, presumably tidally-locked 
gaseous giant planets orbiting only a few stellar radii away 
from their Sun-like host star. By virtue of their slow rota- 
tion (synchronous with their orbital periods, <~ a few days), 
high atmospheric temperatures (r > 1000 K) and permanent 
day-side hemispheric forcing, hot Jupiters are laboratories for 
the study of atmospheric dynamics in a regime that is ab- 
sent from the Solar System (see Showman et al. 2008 and 
Showman et al. 2010 for reviews). In recent years, consid- 
erable progress has been made in observationally characteriz- 
ing the atmospheres of hot Jupiters via a combination of sec- 
ondary eclipse, transmission spectrum and phase curve mea- 
surements (see Deming 2008 and Charbonneau 2009 for re- 
views). In parallel with this observational progress, theoreti- 
cal modeling of hot Jupiter atmospheres has expanded greatly, 
in an attempt to provide robust interpretations of the growing 
data set (see, e.g., Burrows & Orton 2010, Showman et al. 
2010 and Baraffe et al. 2010 for reviews). 

One of the main interests in studying hot Jupiter atmo- 
spheres lies in understanding their thermal and dynamical 
responses to the unusual forcing conditions they are experi- 
encing, with an atmospheric circulation pattern that is likely 
different from anything known in the Solar System. How- 
ever, with this new regime also comes the possibility that new 
physics is at play in these extreme atmospheres (e.g., Menou 
& Rauscher 2010). In this work, we investigate the possibility 
that magnetic drag on atmospheric motions provides an effec- 
tive frictional mechanism limiting the asymptotic speeds of 
winds in hot Jupiter atmospheres. This is particularly impor- 
tant as the fast (transonic) speeds reached by winds in a vari- 
ety of drag-free atmospheric models for this class of planets 
(Dobbs-Dixon & Lin 2008, Dobbs-Dixon et al. 2010; Show- 
man et al. 2009; Rauscher & Menou 2010) raise issues about 
compressibility, shocks and associated energy conservation 
for the models (Goodman 2009; Rauscher & Menou 2010). 
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Atmospheric motions are driven by pressure-gradient 
forces arising from differential heating of the atmosphere. A 
small fraction of the atmospheric "available" enthalpy is con- 
tinuously converted into kinetic energy of the atmospheric 
motions, which is itself continuously dissipated by friction 4 
(Lorenz 1955, Pearce 1978, Marquet 1991, Goodman 2009). 
In steady-state, asymptotic wind speeds are thus reached 
through a detailed balance between continuous thermal forc- 
ing and sustained friction. While the source of wind friction 
on the Earth, and other Solar System terrestrial planets by ex- 
tension, is understood to be largely associated with surface 
drag, the origin of friction in the atmospheres of gaseous gi- 
ant planets remains a major open question in atmospheric sci- 
ence, even in the Solar System (e.g., Schneider & Liu 2009; 
Liu et al. 2008; Showman et al. 2010). Identifying domi- 
nant sources of internal friction in gaseous giant planet atmo- 
spheres can thus be as important as adequately modeling their 
sources of thermal forcing. 

We show here that magnetic drag on weakly-ionized winds 
in the predominantly neutral atmospheres of hot Jupiters, 
which arises from wind interaction with the magnetic field 
generated in the planet's bulk interior, is a plausible source 
of sizable friction, that may need to be accounted for in at- 
mospheric circulation models of hot Jupiters. We note that, 
while this manuscript was being prepared, Batygin & Steven- 
son (2010) completed a study of the closely related ohmic 
dissipation process associated with the currents induced by 
magnetic drag. Here, we focus on the role of magnetic drag 
on atmospheric winds and defer a study of ohmic dissipation 
to future work. 

2. ATMOSPHERIC CURRENTS AND ION DRAG 

2.1. Model Atmosphere 

Our study is based on the three dimensional hot Jupiter at- 
mospheric circulation model of HD 209458b presented by 
Rauscher & Menou (2010). The model provides, at each 
point in the 3D atmosphere, values for the pressure, tem- 
perature, zonal (east-west) and meridional (north-south) wind 
velocity. Throughout this analysis, we use the model atmo- 

4 In this context, friction typically refers to a dissipative process that is 
much more efficient than the microscopic viscosity of the atmospheric gas. 
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sphere obtained after 1450 planetary days of integration. The 
corresponding atmospheric structure and wind pattern are de- 
scribed in detail in Rauscher & Menou (2010). 

The model extends vertically from 1 mbar at the top to 220 
bar at the bottom. Temperatures are typically around 1800 K 
in much of the deep atmosphere, while they vary noticeably 
between day and night near the model top, with lows of about 
500 K and highs of about 1500 K. From the pressure and tem- 
perature, we infer local gas densities using the ideal gas law, 
p = pm^ p/kT ', with mean molecular weight p = 2.33m p . Den- 
sities in the atmosphere range from about 10~ 3 g cm" 3 at the 
model bottom to ~ 10~ 8 - 10~ 7 g cm -3 at the top. 

2.2. Ionization 

High in the atmosphere (say, at nanobar levels), UV photo- 
ionization is important in determining the atmospheric ioniza- 
tion level (e.g., Murray-Clay et al. 2009), but in the relatively 
dense levels modeled here, thermal ionization is expected to 
dominate the ionization balance. 5 

At the temperatures of interest in our model atmosphere, the 
main source of free electrons is provided by thermal ioniza- 
tion of alkali metals with low first-ionization potentials: Na, 
Al, and K. Under these conditions, the mean ion mass is on 
the order of m, rj 30m p (see e.g. Draine et al. 1983). For 
simplicity, we choose to approximate Saha's equation for ion- 
ization balance with a formulation that only accounts for the 
ionization of potassium (Balbus & Hawley 2000), 
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where n e and n„ are the number densities of electrons and 
of neutrals, respectively (in cm" 3 ), a K is the potassium abun- 
dance, and T is the temperature in K. Equation ([T} is a valid 
approximation only as long the resulting ionization fraction, 
x e , remains much smaller than the abundance of potassium, 
aK- We have verified that this condition is reasonably well 
satisfied for the atmospheric conditions of interest here, with 
x e reaching at most ~ 10" 9 in a few localized regions and 
taking much smaller values (x e ~ 10~ 10 -10~ 14 ) in the rest of 
the atmosphere. Throughout our analysis, we assume a near 
solar abundance of potassium, a^ = 10~ 7 . This is not a crit- 
ical model assumption given the much more important, ex- 
ponential dependence of x e with temperature. We also esti- 
mate that, even at T ~ 1800 K, the free electron contribution 
from sodium, which is ~ 17 times more abundant than potas- 
sium at solar composition, approaches only marginally that 
of potassium. While a more complete Saha equation solution 
would be a clear improvement upon the very simple approach 
adopted here, it would not qualitatively alter our main conclu- 
sions about the role of magnetic drag. 

We assume conditions of gas neutrality, i.e. n e = n,-, where 
«, is the ionic number density. The coupling between charged 
(i, e) and neutral particles depends on their rate of collisions, 
(erv) ( e , for which we adopt the expressions from Draine et al. 
(1983) 
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X-rays could provide a source of additional non-thermal electrons at and 
above mbar pressure levels, which is ignored from our calculation (J. Good- 
man, priv. communication). 
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The electrical conductivity, <j e , and the corresponding resis- 
tivity, r/, are respectively given by 
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2.3. Non-ideal MHD regime 

Before proceeding with a calculation of induced currents, 
we need to evaluate the relative importance of various non- 
ideal MHD effects. We begin by considering the full induc- 
tion equation 
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where the last three terms are the ohmic, Hall and ambipo- 
lar diffusion terms, respectively. In the above equation, the 
current density 

J=^(VxB), (5) 

47T 

v is the velocity of the neutrals, p, is the ion mass density and 
7 is the drag coefficient, which is given by 



7 : 



nii + p 



(6) 



We estimate the relative importance of the three non-ideal 
MHD terms by selecting relevant velocity and length scales 
for the atmospheric flow. For a representative induction ve- 
locity, we take the sound speed c s , since it is representative 
of the wind speeds achieved in the modeled transonic flow. 
Following Rauscher & Menou (2009), we use the adiabatic 
value c s = y/{TkT / pm p ), with an adiabatic index V = 1.47. 
For now, we adopt a pressure scale height, H p , as a represen- 
tative length scale for our comparison of the non-ideal MHD 
terms, even though we argue below that using the resistive 
scale height may be a more appropriate choice. With these 
choices, the magnetic Reynolds number can be evaluated as 
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where R' = R/ p = 3.78 x 10 7 erg/g/K is the effective gas 
constant and g = 9.42 m/s 2 is the gravitational acceleration of 
the planet. 

We evaluate the relative magnitude of the Ohmic (Ohm), 
Hall (Hall) and Ambipolar diffusion (Amb) terms with respect 
to the Induction term (Ind), using the same notation as Balbus 
& Terquem (2001; see also Sano & Stone 2002) for conve- 
nience, 
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with a velocity scale, v = c s , and length scale, L = H p . 
We computed these ratios throughout the modeled atmo- 
sphere, assuming a nominal field strength Bq = 3 G, and found 
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FIG. 1. — Cylindrical maps of the magnetic Reynolds number, Logio(R,„), on two different pressure levels in the model atmosphere. The sub-stellar point is 
centered at longitude and latitude zero. Left: highest model level (p = 1 mbar); Right: moderately deep level (p = 2 bar). While very large day-night variations in 
R m occur high in the atmosphere, R m values are much more uniform deeper in the atmosphere. At all levels, R m only marginally exceeds unity. 



that the Ohmic term is typically larger than the Hall term, by 
about 12 orders of magnitude in the deepest levels, down to 
about 3 orders of magnitude in the highest level modeled. The 
Ohmic term is also found to dominate over the Ambipolar dif- 
fusion term, although to a lesser extent. It is about 5 orders 
of magnitude larger in the deepest levels, down to a factor of 
only a few in the uppermost modeled level. We conclude that 
the Hall term can be safely neglected throughout the mod- 
eled atmosphere, while the Ambipolar diffusion term could in 
principle become important in the more tenuous high atmo- 
sphere, at and above mbar pressure levels, especially if mag- 
netic field strengths B > 10 G are considered. Based on these 
results and in the interest of simplicity, we focus the rest of 
our analysis on the purely resistive MHD regime. 

Figure 1 shows cylindrical maps of the magnetic Reynolds 
number on two different pressure levels in our model atmo- 
sphere: the highest level, at p = 1 mbar (left panel), and a 
deeper level with p ~ 2 bar (right). These magnetic Reynolds 
number maps largely trace the variations in conductivity at a 
given pressure level. High in the atmosphere, despite strong 
winds, the large temperature difference between the day and 
the night side leads to many orders of magnitude variations in 
R m values. By contrast, deeper in, temperature variations on 
a given pressure level are modest and so are the variations in 

Rrn ■ 

For the treatment of magnetic drag presented in j p.4| to be 
valid, it is important that R m <C 1 in the model atmosphere. 
Indeed, an atmospheric flow with R m > 1 is coupled to the 
magnetic field well enough that, in principle, it could gen- 
erate its own magnetic field via dynamo action, leading to a 
situation that could significantly complicate any discussion of 
magnetic drag. Examining values of R„, throughout our model 
atmosphere, we find that they can exceed unity in particular 
regions of the atmosphere (as shown, e.g., in Figure 1), but 
only marginally so. 

However, it is also likely that the simple definition of R„, 
adopted in Eq. (|7|i effectively overestimates the ability of 
the model atmosphere to generate its own dynamo. Indeed, 
the length scale for dimensional analysis in a geometrically 
anisotropic system like an atmosphere must be chosen care- 
fully. Examination of the specific toroidal component of the 
induction equation used below (see Eq. |9)) shows that ver- 
tical gradients generally dominate over latitudinal ones. The 
choice of the vertical gradient scale itself is subtle since it is 
different for pressure (H p ) and resistivity (H rj ), with a gener- 



ally shorter resistivity scale from its exponential dependence 
on temperature. Based on the leading order vertical gradient 
terms entering Eq. d9]l, this suggests that H n may be the rele- 
vant scale to define R,„ in our problem, rather than H p . Indeed, 
H v is generally smaller than H p in our model atmosphere, by 
a factor of a few typically on the dayside and a factor ~ 10 to 
100 on the nightside in the upper levels, and by a factor of a 
few deeper in. Based on these arguments, we assume for the 
remainder of this analysis that there is no self-induced atmo- 
spheric dynamo, so that the planetary magnetic field, originat- 
ing from deep currents in the bulk interior, is the only relevant 
one. 

2.4. Induced currents and magnetic drag 

Since our goal in this work is to establish whether or not the 
typical magnitude of magnetic drag in hot Jupiter atmospheres 
is sufficient to appreciably affect their circulation patterns, we 
adopt the simplest approach possible and make a number of 
simplifying assumptions. We assume the planetary magnetic 
field to have a dipolar geometry, with a dipole axis coincident 
with the rotation axis. The surface magnetic field strengths 
expected for hot Jupiters are rather uncertain. Jupiter has a 
field strength varying from about 4 G at the equator to about 
15 G at the poles. In our calculations, we assume a fiducial 
surface field strength Bq = 3 G (at the magnetic pole) and dis- 
cuss the effects of larger field values whenever appropriate. 

We adopt a formalism that is closely related to that de- 
scribed by Liu et al. (2008). In particular, we strictly focus 
our analysis on the zonal (= east-west) component of the at- 
mospheric flow and the drag experienced by this flow as it 
generates a toroidal field component, and associated currents, 
from the purely poloidal planetary magnetic field. The gener- 
ation of the toroidal component of the magnetic field is gov- 
erned by the resistive induction equation 
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where Q = v^r" 1 sin" 9 in spherical coordinates (r, 9, (f>). We 
seek a steady-state solution to this induction equation on the 
assumption that it represents a reasonable value for the level 
of drag expected in our model atmosphere with prescribed 
wind speeds. For simplicity, we entirely neglect the poloidal 
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component of the induction equation, which would lead to 
magnetic drag on the meridional (=north-south) component 
of the circulation. While this assumption may be justifiable in 
much of the model atmosphere, where the dominant circula- 
tion is zonal, it is likely a poor approximation high up in the 
atmosphere, where the circulation pattern is away from the 
substellar point and thus equally meridional as zonal. 

In our three-dimensional model atmosphere, resistivity 
varies strongly in the vertical but it also does so hor- 
izontally, in the upper atmosphere (see, e.g., the left 
panel of Fig. 1). It is thus not a priori clear whether 
the second-to-last term dominates over the last term in 
Eq. The relative magnitude of these two terms de- 

pends on the relative magnitude of the vertical gradient scale 
Hrj(r) = \rj(r,8,(j))/(dr)(r,6,(f>)/dr)\g ^ and the horizontal gra- 
dient scale H„(0) = (1 /r) \r](r,0, $)j{difyr, 0, <WA^)U for re- 
sistivity. We computed the ratio H„(f) / HJff) in our model 
atmosphere and found it to be < 0.1 over the vast majority of 
the atmospheric domain. As a result, the last term in Eq. |9]l 
is small and can be dropped from the analysis. The domain 
in which horizontal resistivity gradients can approach vertical 
gradient values is about 1 % of the total, and about 2/3 of this 
domain resides in the model uppermost layers, at p < 10 mbar. 
For these upper levels, our estimates for the magnetic drag 
may thus be inaccurate in some regions because the strong 
horizontal resistivity gradients can no longer be ignored. 

Under the various assumptions made so far, the meridional 
component of the current induced by the zonal flow is domi- 



nant and can be evaluated as (Liu et al 
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where the last term is associated with a boundary current in 
the uppermost modeled level, jg(R,8,(f>). Eq. (JXTJ expresses 
that the local current scales as an integral of the zonal flow 
above the level of interest. Lacking information about the 
nature of currents possibly flowing from regions above the 
modeled atmospheric layers, we set this boundary current to 
zero for simplicity. This represents an important source of un- 
certainty in our modeling approach, but we note that, unless 
near cancellations occur, additional boundary currents could 
in principle contribute to even stronger magnetic drag than 
estimated here. 

We evaluate the magnitude of the magnetic drag on the 
(mostly neutral) flow by using a standard formulation for "ion 
drag" based on the bulk Lorentz force experienced by the 
ionic component (see, e.g., Zhu et al. 2005), 



dy 1 

p— oc - j xB , 



(12) 



dt c 

where j is the current induced in the atmosphere by the zonal 
flow. From this, we deduce the typical zonal drag time over 
which the zonal flow would be brought to a halt, 
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in the absence of any other forces. 

An order of magnitude estimate for the drag can be sim- 
ply derived by making the approximation jg ~ cv $B /(Anrj), 
which yields 
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FIG. 2. — Cylindrical maps of the zonal drag time, Logioiy<irag (s)], on four 
different pressure levels in the model atmosphere. The sub-stellar point is 
centered at longitude and latitude zero. From top to bottom: highest model 
level, at p = 1 mbar; p = 40 mbar level; p = 0.5 bar level; and last a moder- 
ately deep level, at p = 2 bar. High in the atmosphere (top two panels), drag 
times vary very strongly between the day and the night sides. Drag times 
also become very long in the vicinity of the equator, by virtue of the aligned 
dipolar field geometry adopted. 
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where B 3 = 5/(3 G), p_ 5 
?7/(10 13 cm 2 s _1 ) are typical values, and 7r - 
gle between jg and B. 

Fig. 2 displays maps of the magnetic drag time from the 
detailed modeling of Eq. (13 1, at four depths within the at- 



mosphere. At deeper levels, the drag time is relatively long, 
Tdrag > 1 7 — 1 8 sec, with only a few localized pockets with 
shorter drag times, corresponding to regions of the flow with 
low velocities. Higher in the atmosphere, the magnetic drag 
time spans a very large range of values, with very long drag 
times on the nightside and considerably shorter drag times on 
the dayside. While magnetic drag is probably negligible on 
the night side, the short drag times ~ 10 4 - 10 6 sec on the 
day side suggest that magnetic drag could be dynamically im- 
portant in these regions. In addition to the day-night asym- 
metry in drag times high up in the atmosphere, Fig. 2 shows 
that zonal drag times are consistently long along the equator. 
This geometrical effect, which results from the aligned dipo- 
lar field geometry assumed in our analysis (leading to B r = at 
the equator), could have interesting consequences for the at- 
mospheric flow, such as promoting barotropic shear instabil- 
ities in the equatorial regions (see, e.g., Menou & Rauscher 
2009). This geometrical effect may persist for dipolar field 
geometries with a small amount of misalignment. 

2.5. Preliminary models with Rayleigh drag 

Since the shortest drag times shown in Fig. 2 are compa- 
rable to a representative zonal advection time in our model 
atmosphere, r a( j v ~ R P /c s ~5 x 10 4 s, one may expect such 
levels of magnetic drag to affect the atmospheric flow non- 
trivially, even in the presence of continuous thermal forcing 
of the circulation. 

To test this hypothesis, we have run modified versions of 
the drag-free model described in Rauscher & Menou (2010) 
in which a vertically-varying but horizontally-uniform level 
of Rayleigh drag (oc -v/T,j rag ; see Menou & Rauscher 2009) 
is applied throughout the modeled flow. For simplicity, an 
equal amount of drag is applied to the zonal and the merid- 
ional velocity fields and only average, representative values of 
Tdrag for each pressure level are employed. Assuming a field 
strength Bq = 3 G, we evaluated representative horizontally- 
averaged values of Td rag ~ 6 x 10 6 s at 1 mbar (model top), 
~ 4 x 10 7 s at 40 mbar, ~ 7 x 10 7 s at 2 bar and ~ 8 x 10 8 s 
at 220 bar (model bottom), with a logarithmic dependence of 
Tdrag on pressure in between. 

We have run models for the nominal field strength Bo = 3 G 
as well as models with drag times uniformly reduced by a fac- 
tor 10 and 100 (corresponding to field strengths B ~ 3 x Bq 
and 10 x Bo)- Each model was run for 500 planet days ^rota- 
tion periods). In the model with B = Bq, a very similar atmo- 
spheric flow structure to that described in Rauscher & Menou 
(2010) is obtained, with zonal wind speeds slightly reduced 
at the ~ 5-10% level (Fig. [3] top panel). In the model with 
B = 3x Bq, the atmospheric flow is again qualitatively similar 
but the superrotating equatorial jet is shallower and confined 
to levels above ~ 1 bar, while zonal wind speeds are typically 
reduced by a factor ~ 30% (Fig. [3] middle panel). Finally, in 
the model with B = 10 x Bq, we witness significant changes in 
the structure of the atmospheric flow. Zonally-averaged zonal 
wind speeds in the equatorial jet are reduced by a factor ~ 3 
while stronger counter jets dominate at mid-latitudes (Fig. [3] 
bottom panel). Nevertheless, peak velocities in the upper at- 
mosphere remain supersonic. While these three simple runs 



are only illustrative, since they do not differentiate zonal and 
meridional drag or include the strong horizontal dependence 
of the drag (Fig. 2), they clearly suggest that magnetic drag of 
the typical magnitude evaluated in this work can have a no- 
ticeable, and possibly important, impact on the atmospheric 
flow of hot Jupiters. 

3. SUMMARY AND DISCUSSION 

Circulation models for hot Jupiter atmospheres have gained 
substantial interest in recent years thanks to the direct detec- 
tions of such atmospheres, which indicate a possible role for 
winds in shaping the emergent properties of these planets. All 
existing calculations have been purely hydrodynamical in na- 
ture, neglecting the possibility that interactions between the 
flow and the planetary magnetic field could influence the cir- 
culation. However, given the high atmospheric temperatures 
of these planets, the small fraction of the fluid that is ionized 
may be sufficiently large to effectively couple the mostly neu- 
tral flow to the planetary magnetic field. 

In this work, we have evaluated the typical magnitude of 
the magnetic drag exerted on hot Jupiter atmospheric winds, 
using the flow pattern obtained in the 3D circulation model of 
Rauscher & Menou (2010). Our results indicate that typical 
magnetic drag stopping times are a strong function of depth, 
longitude and latitude. In the upper model atmosphere, tem- 
perature differences between the day and night side cause ex- 
treme variations of the magnitude of the drag, which could 
be dynamically important only on the day side. At all levels, 
zonal drag also varies strongly with latitude in the vicinity 
of the equator, by virtue of the aligned dipole field geometry 
adopted. 

By adopting a simplified treatment of ionization balance, 
focusing on the zonal component of drag only, ignoring up- 
per boundary currents or the role of ambipolar diffusion, and 
neglecting horizontal variations in resistivity, the calculations 
presented in this work are clearly of limited scope. Neverthe- 
less, they appear sufficient to make the case that dynamically 
interesting levels of magnetic drag may have to be accounted 
for in circulation models for hot Jupiter atmospheres. This 
conclusion is supported by a simple comparison between a 
drag-free circulation model and a few additional models with 
Rayleigh drag applied at levels commensurate with that ex- 
pected from magnetic drag, which show dynamically interest- 
ing consequences for the flow, especially for large magnetic 
field values. 

A better assessment of the role of magnetic drag in hot 
Jupiter atmospheres may require iterated solutions to the 
multi-dimensional induction equation, together with Saha's 
equation, as the time-dependent hydrodynamical flow is 
evolved. Such work will be important to re-evaluate the role 
of deep ohmic dissipation in possibly inflating hot Jupiters 
(Batygin & Stevenson 2010) and to elucidate the plausible di- 
versity of atmospheric behaviors expected from a variety of 
field strengths and geometries. The strong dependence of re- 
sistivity on atmospheric temperature could, in principle, lead 
to different classes of atmospheric behaviors as a function of 
mean orbital separation, and it will likely lead to an inter- 
esting phenomenology for eccentric close-in planets, as they 
experience time-variable levels of magnetic drag (and ohmic 
dissipation). 

We thank Adam Burrows, Jeremy Goodman and Peter Gol- 
dreich for useful discussions and encouragements. 
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FIG. 3. — Zonal average of the zonal wind in m s , as a function of latitude and depth, in two atmospheric circulation models with vertically-varying and 
horizontally-uniform levels of Rayleigh drag. Yellow lines separate regions of positive (eastward) flow and negative (westward) flow. These contour plots are 
meant to be compared to the equivalent version for the drag-free model described in Rauscher & Menou (2010). Top: Rayleigh drag times calibrated to match the 
typical magnetic drag times expected for a field strength fi~ 3 G; Middle: Rayleigh drag times calibrated to match a magnetic field strength B ~ 10 G. Bottom: 
Rayleigh drag times calibrated to match a magnetic field strength B ~ 30 G. Significant modifications to the atmospheric flow emerge, relative to the drag-free 
case, particularly in terms of the zonal-mean wind speeds (see text for details). 
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